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Abstract 


An algorithm for the optimum design of an internal flow component to 
obtain the maximum pressure rise is presented. Maximum pressure rise in a 
duct with simultaneous turning and diffusion is shown to be related the 
control of flow separation on the passage walls. Such a flow is usually 
associated with downstream conditions that are desirable in turbomachinery 
and propulsion applications to ensure low loss and stable performance. The 
algorithm requires the solution of an "adjoint" problem in addition to the 
"direct" equations governing the flow in a body, which in the present 
analysis are assumed to be the laminar Navier-Stokes equations. Earlier 
studies have usually addressed such problems for the case of inviscid 
and/or irrotational flow. These assumptions may not be valid in flows that 
undergo sharp turning resulting in strong secondary flows and possibly 
separating and recirculating regions. The theoretical framework and 
computational algorithms presented in this study are for the steady Navier- 
Stokes equations. 

A novel procedure is developed for the numerical solution of the 
adjoint equations. This procedure is coupled with a direct solver in a 
design iteration loop, that provides a new shape with a higher pressure 
rise. This procedure is first validated for the design of optimum plane 
diffusers in two-dimensional flow. The direct Navier-Stokes and the 
"adjoint" equations are solved using a finite volume formulation for 
spatial discretization in an artificial compressibility framework. The 
discretized equations are integrated using explicit Runge-Kutta time steps 
to obtain steady-state solutions. It is found that the computational work 
required to solve the "adjoint" problem is of the same order as that 
required to solve the direct problem. It is also found that the procedure 
converges within about ten iterations, and in addition, the number of 
design iterations are not sensitive to the grid used for the calculations. 
This is a significant computational advantage over heuristic design 
procedures based on point by point sensitivity analysis where the work 
increases with the refinement of the grid. 

A simplified version of the above approach is then utilized to design 
ninety degree diffusing bends. The bend inlet is square with intermediate 
and exit cross-sections constrained to be rectangular. The location of bend 
walls is then determined in order to obtain the maximum pressure rise 
through the bend. Calculations were carried out for a mean radius ratio at 
inlet of 2.5 and Reynolds numbers varying from 10f) to 500. While at this 
stage laminar flow is assumed it is shown that a similar approach can be 
conceived for turbulent flow’s. 
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MOTIVATION 


# How to shape internal flow passages 

# Combined turning and diffusing flow 

# Maximize pressure rise 

# Can not assume inviscid or 2-D flow 

OBJECTIVES 

# Develop theoretical framework — laminar 3-D 

# Develop Navier-Stokes and Adjoint solvers 

# Validate Navier-Stokes solver 

(Laminar 90 degree bend, Taylor et al. 1982) 

# Validate Optimization Approach on 
2-D straight diffusers 

# Apply to the design of ninety degree bends 
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'U'i,! — 0 

UjUiJ = -p*i + UUiJj , p* =p/p 

# No slip BC on Tm 

# Dirichlet BC for at Tj and To 



p*UiUids + 
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p*Uin{ds 





J (r*M e ) ~ J(Xm) — + 0(e 2 ) 



# Our “Adjoint” equation: 


V z i, jj Zj t i) r ,i 0 


in fi 


m = - & - ui ) 


q = | (r - P* + (l/2) u i - 2«ywy) 


# Pironneau’s “Adjoint” equation: 

Wii = 0 


vw i} jj + UjWi'j + WjUj'i - q t i = —UjUij 


U); = 0 


in Q 
in fi 
on r 


8J 


= V [ p{s) ( 


du{\ ( dui dwi 


1 l 1 - 2 - 1 ds 

v M \dnj\dn dn 


dw 

dt 


= w (const) + • • • 
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NAVIER-STOKES SOLVER 


p t = -P 2 u it i 
'U'ij “h UjUij — 

a) Artificial Compressibility 

b) Runge-Kutta Time Integration 

c) Finite Volume Discretization 

d) Artificial Dissipation 

e) Local Time Stepping 

e) Implicit Residual Smoothing 
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Figure 4: Geometry of a circular bend with square cross section. 




Figure 5: Streamwise velocities in a bend: a) 0 = 30 degrees, b) © — 60 degrees, c) 0 77.5 degrees. 
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ADJOINT EQUATION SOLVER 



MESH GENERATION 
-#■ Thompson et al. 

&££ %TJT] 0 

vet + y V v = 0 
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Boundary Conditions 

# No-slip be on the walls 

# Zero normal derivatives at exit 

# Streamwise velocity component is specified at 
entrance 

# Zero normal derivatives for remaining velocities 
at entrance 

# Typical bc’s at symettry planes 

=#= Zero second derivatives for pressure (a compu- 
tational be) 



Diffuser Profile History 





Wall Shear Stress 


Skin Friction History 



Re=200, 61 by 31 grid 

o : N=1 

n : N=2 
A : N=5 
* : N=10 
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Pressure Rise 


Pressure Rise History 



Re=200, 61 by 31 grid 

O : Area-averaged pressure rise 
A : Flow-averaged pressure rise 
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Pres 


Actual Pressure Rise 


Ideal Pressure Rise 



Reynolds Number 


O • Optimal diffusers 
A : Straight diverging diffusers 


2 






Grid Study 





T 

2 


Re=200 

O : 31 by 16 
n : 61 by 31 
* : 121 by 61 




The error in the location of the optimal diffuser 
profile corresponding to a 2 percent error in the 
total pressure rise. The optimum shape lies between 
the high and low y-values shown in the graph 
(Re=500, grid size is 61 by 21, and L/Wi = 3). 
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Sketch of the inlet header 



Figure 3: A representative section of a three-dimensional diffuser. Flow enters at upstream boundary T/ 
and exits at downstream boundary To- The walls to be shaped are Im- 
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ISSUES 


# Geometry Constraints 

In general: Move walls by ep(s) along the 

normal direction, everywhere 

New shape may not satisfy 

— specified mean passage location 

— specified cross-sectional shape 

— overall system geometry 

# Present Work 

- No correction on side walls (z = 0, z = 2r max ) 

- Apply mid-plane (z = z m ax /2) correction 
to all £ locations 

- Hence all cross-sections are rectangular 

# Laminar flow results (Re < 500) 
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Governing Equations: 


u j u ijj = ”"P,t "b l/u i,jj 


Design Objective: 

Maximize Static Pressure Rise 



Cabuk and Modi, 1990 




= 6 
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DESIGN ALGORITHM 


1: Choose an initial shape. 

Generate the computational grid. 

3: Solve the N-S equations. 

4: Compute shear stress on the walls. 

5: Compare wall shear stress to target dis- 
tribution and determine the amount of 
boundary movement p(s). 

6: Update the shape. 

7: Go to step (2) 
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Normalized wall 
shear stress 



Wall shear stress along the outer wall, 
Re=100 

O • Optimum diffusing bend 
A : Elliptic diffusing bend 
Dashed Curve : Target distribution 
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Normalized wall 
shear stress 


1.4 



Arclength along the bend 


Wall shear stress along the inner wall, 
Re=100 

O • Optimum diffusing bend 
A : Elliptic diffusing bend 
Dashed Curve : Target distribution 
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Optimum Shapes 
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Performance = f(shape) 

^ Possible Applications: 

- 90 Degree Bend 

- Turn Around Ducts 
— Transition Ducts 

- S-shaped Ducts 

— Straight or Curved Diffusers 
— Turbine Blades 

- Engine Inlets 

- Turning Vanes 
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CONCLUSIONS 


Theory 

# Theoretical Framework for Design 
with Navier-Stokes equations 

# Determine p(s) from Direct-f- Adjoint 
or from Direct alone 

Computational 

# Direct and Adjoint Solvers Validated 
for Plane Diffusers 

# Design of 90 degree Bend with Specified 
Cross-section, Max. Ap 

# Number of Design Cycles < 10 

# “Flow” Interpretation of Adjoint Problem 

Future Plan 

# Apply to 3-D turbulent flow 

# Specify mean line, vary cross-section 

# Other objectives : Min. Distortion 
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